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Ov 

• The collapse of a massless scalar field in the Brans-Dicke theory of gravitation is studied in the 

analysis of both analytical solution and numerical one. By conformally transforming the Roberts's 
solution into the Brans-Dicke frame, we find for u > —3/2 that a continuous self-similarity continues 
and that the critical exponent does depend on u). By conformally transforming the Choptuik's 
solution into the Brans-Dicke frame, we find for uj > —3/2 that at the critical solution shows discrete 
qq . self-similarity, however, the critical exponent depends strongly on uj while the echoing parameter 

' weakly on it. 

7— I . 

> 1 I. INTRODUCTION 

\o : 

IT) , 

Critical phenomena in black hole formation found by Choptuik |1| has renewed interest in the classical general 

£f) • relativistic black hole formation. Among the most interesting things in Choptuik's results are 

(1) Power law: Black hole mass exhibits a power law 

\o : 

M BH ^{p~p*f, (1.1) 

where p is a parameter which characterizes the strength of initial condition and p* is the threshold value. 

(2) Echoing: Configurations sufficiently close to critical show a discrete homotheticity (or scale invariance) 

0( P -A,t-A)^(p,t), (1.2) 

where p and r are logarithms of proper radius r and central proper time t. Echoing means that the features of critical 
rS | solution are repeated on ever decreasing time length scales. 

(3) Universality: The exponent (5 ~ 0.37 and echoing parameter A are independent of any choice of initial data p. 

It should be noted, however, that the universality mentioned there refers to the independence from the initial 
condition of the matter field considered. To investigate the dependence of the model of the matter, the collapse of 
gravitational wave J|] and radiation fluid collapse || were examined. For vacuum gravity, [3 ~ 0.37 and A ~ 0.6, 
while for radiation fluid [3 ~ 0.36 with A being arbitrary (i.e. continuous self-similarity). These results excite the 
expectation that the critical exponent (3 may be universal among matter fields having massless property, although 
quantum effects will destroy the phenomena against expectation [EJ. 

In this paper, we will investigate another model dependence: dependence on the theory of gravitation. To see this we 
take the Brans-Dicke theory of gravitation |q] for its simplicity. The theory contains a parameter uj which controls the 
strength of gravity, and in the limit uj — > oo the theory reduces to the Einstein theory. We consider the Brans-Dicke 
theory in vacuum. The scalar field examined here is thus the Brans-Dicke scalar field. 

The collapse of a massless scalar field in the Brans-Dicke theory is also interesting in the light of structural stability 
of Choptuik's solution, which is our next problem. 
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Our paper is organized as follows. In sec. 2, Roberts's self-similar solution in the Brans-Dicke theory is studied 
to introduce the notation and the detail of the conformal transformation. In sec. 3, basic equations and numerical 
procedures are given. We reproduce Choptuik solution. In sec. 4, numerical results for Brans-Dicke theory are given. 
Sec. 5 is devoted to summary. 



II. ROBERTS'S SOLUTION IN BRANS-DICKE THEORY 

Roberts || found solutions of the Einstein equations which describe continuous self-similar spherical collapse of 
a massless scalar field and he studied these solutions in the context of counter-examples to the cosmic censorship. 
Subsequntly, the solutions were studied as an analytical example of critical behavior. We shall try to find the 
counterparts in the Brans-Dicke theory of Roberts's solutions by conformal transformation. 

The action of Brans-Dicke theory without matter term is 



By the conformal transformation 



S= I AfjtM-^Vf) 2 ). (2.1) 



9ab = $ 1 9ab, 



$ = exp(0/V2w + 3), (2.2) 



the action Eq.(2.1) is reduced to that of a massless scalar field coupled to the Einstein gravity where the solutions are 
known: 

#Xy/=g{R - i(V^) 2 ). (2.3) 



Therefore we simply replace the Roberts's solutions by Eq.(2.2) to find the counterparts in the Brans-Dicke theory. 
The Roberts's solutions are written as 

ds 2 = -dudv + ~(u- (l-p)v)(u - (l+p)v)dn 2 , (2.4) 

where p > is an arbitrary constant. Then those in the Brans-Dicke theory are given by 

dS 2 = ^ ) ±a dudv + ~(u- (l~p)v) 1±a (u- {l+p)v) x * a d&, (2.6) 



<S>=( U - { '- p \ V Y a , (2.7) 
\u-(l+p)vJ y ' 



where a = \ j \J2uj + 3. Since the Ricci scalar takes the form 

£ = $.R + 3V 2 $-^(V$) 2 , (2.8) 

and R — 8[(u — (1 — p)v))(u — (1 + p)v)]~ 2 p 2 uv, we find that a singularity at u = (1 — p)v for v > and u = (1 +p)v 
for v < 0. In the following we shall assume 1 ± a > in order for the singularity to locate at the center. As in the 
case of the Einstein gravity, structure of the singularity depends on p. It is easily found that for p > 1 it is timelike 
in the region v < and spacelike v > 0, null for p — 1, timelike for p < 1. 

By calculating the change of area with respect to derivative along advanced null coordinate, it is found that an 
apparent horizon is located at 

1 -v 2 

u = —^v, (2.9) 
1 ± otp 

and that for p < 1 no apparent horizon exists. Hence the parameter p corresponds to the control parameter in this 
example of the critical behavior and p = 1 the point of the criticality. 
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Calculation of the black hole mass M is also straightforward. It is done by calculating the square root of the proper 
area of an apparent horizon. After a little algebra it is found that M is given by 

(1 ± q)( W (1 zp a) (l± a )/>(l +p)(lTa)/2 

4(-l±op) (P ' ' 1 j 

Thus near the criticality M behaves as M oc (p — 1)( 1±Q )/ 2 . It is the power law of the critical behavior with the 
critical exponent 

0= - + — , (2.11) 

which clearly depends on a; and certainly reduces to the result of Einstein theory in the limit to — > oo. This result 
strongly indicates that the critical phenomena depend on the parameter lj. Indeed, it will be confirmed in the 
subsequent numerical analysis. 



III. COLLAPSE OF BRANS-DICKE SCALAR FIELD IN DOUBLE-NULL COORDINATES 

Next we study numerically the collapse of the Brans-Dicke scalar field. We treat the problem with the action 



Eq.(2.3) in terms of null initial value formulation. Since in the null initial value formulation the grid points are tied 
to ingoing light rays, overall size of the grid becomes smaller as system evolves and the resolution improves. Adaptive 
mesh refinement algorithm used by Choptuik is not always necessary. 
We take the following line element 

ds 2 = -a(u, vfdudv + r(u, vfdti 2 . (3.1) 

Double-null coordinates are not unique. There remain two degrees of gauge freedom which correspond to redefining 
u and v. We fix one of them by requiring u = v on axis r = 0. The remaining gauge freedom will be fixed by the 
initial condition. 

The equations of motion that derived from Eq.( |2.3| ) are 

1 2 , . 

rr uv + r u r v + -a = 0, (3.2) 

a~ 1 a uv - a~ 2 a u a v - r~ 2 r u r v - -^r~ 2 a 2 + ^> u <t>v = 0, (3.3) 

r^uv + r u 4>v + r v 4> u = 0. (3.4) 

And constraint equations are 

r uu - 2aT x a u r u + ^r<p 2 u = 0, (3.5) 

r vv - 2a~ 1 a v r v + ^r<p 2 v = 0, (3.6) 

where r u = dr/du, for example. These equations are solved numerically. Since the usual iteration scheme for solving 
these equations does not work because of their nonlinearity, we adopt the first order scheme developed by Hamade 
and Stewart 0. 

We introduce the following new variables 

d = — ,/ = r u ,g = r v ,p = 4> u ,q = <j> v . (3.7) 
a 



The system of equations (3.2-3.6]) can then be converted to a first order system: 



U+r- 1 (fg + 


\a 2 ) = 0, 


(3.8) 


- 2 (/5 + J« 2 H 


- \pq = 0, 


(3.9) 


Pv +r~ 1 (fq 


+ flp) = o, 


(3.10) 


In +r~ 1 (fq 


+ 9P) = 0, 


(3.11) 


g v - 2dg + 


\rq 2 = 0, 


(3.12) 
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and the equations denning variables 



a v — ad = 0, 
r v - g = 0, 
<f> v - q = 0. 

We have to spe cify the boundary conditions on axis, u = v. 
i.e., / = — g. Eq. p.q ) will be regular on axis if and only if a - 
value for a and 4> are obtained by requiring a r — (j) r — there. 

We need to give initial data d and q on the initial surface u 
second paragraph is the choice of d on u 
(or <t>) on u — 0, we typically take 



(3.13) 
(3.14) 
(3.15) 

= there and this implies r u + r v = 0, 

The boundary 



Of course 

; 2g. Eq.(3.1C) implies p = q on axis 



0. The remaining gauge freedom mentioned in the 
0. This is arbitrary, and we choose d = 0. As for the initial condition of q 



4>{u = 0,v) = 1 + (/> Q v 2 exp(-(« - v a ) 2 /vl). 



(3.16) 



We then integrat e out ward from axis Eq.( |3.13D to obtain a, Eq.( |3.14| ) to obtain r, Eq. fl3.12| ) to obtain g, Eq.(|3~8|) to 
obtain / and Eq.( 3.1C ) to obtain p. Thus we set initial conditions. 

The integration in the it-direction is done using an explicit difference algorithm. The integration in the u-direction 
is done using an implicit algorithm, which is necessary to ensure stability. However, the integration can be made 
explicit. Details of algorithm are given in (?|.| 

In Fig.l, we plot the scalar field at the center with marginally subcritical initial parameter as a function of logarithm 
of the central proper time r = — log(t* — t), where t* is the time when an infinitesimal black hole forms. We find that 
the scalar field oscillates with period A ~ 3.43 in r-coordinate, in agreement with Choptuik. 



Fig. 2 shows the black hole mass for marginally supercritical evolution as a function of <pQ — 0q with 
critical parameter. The least square fit shows the exponent (3 ~ 0.38, in agreement with Choptuik. 
These results confirm the logical consistency and accuracy of our numerical code. 



being the 



IV. CRITICAL PHENOMENA IN THE BRANS-DICKE THEORY OF GRAVITATION 



Now we study the behavior of a scalar field for various values of oj, 1000 > 2uj + 3 > by c onformally transforming 
the Choptuik's solution reproduced in the previous section. This is done by the relation (2.2). The results are shown 
in Figures and Table. Fig. 3 shows the scalar filed at the center as a function of logarithm of the central proper time 
t = — \og(t* — t) for several values of ui. The scalar field oscillates greatly for smaller ui. Fig. 4 shows the black hole 
mass for marginally supercritical evolution as a function of initial parameter for 2lu + 3 = 2. The least square fit gives 
the exponent (3 ~ 0.43. Table summarizes the numerical results. As can be seen the exponent [3 depends strongly 
on oj, hence it is not universal quantity. This is the first numerical evidence of non- universality Q for the collapse of 
single matter content. On the other hand the scaling parameter A depends weakly on uj. 



Table Brans-Dicke parameter, critical exponents, and scaling parameters. 



2w+3 


P 


A 


0.50 


0.50 


3.41 


1.0 


0.46 


3.43 


2.0 


0.43 


3.45 


5.0 


0.41 


3.46 


10 


0.38 


3.46 


50 


0.38 


3.47 


200 


0.38 


3.47 


1000 


0.38 


3.47 



However, slight modification is necessary to ensure that the system has second-order accuracy; Eq.(3.5) in their paper should 



be replaced with z n 



+ 2, 



+ ^(G(y„, z n ) + G(y w , z w ))^j . We thank T.Harada for pointing this out to us. 
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It is to be noted that critical phenomena in the Brans-Dicke theory does emerge in discrete self-similar manner. 
Thus stability analysis around a continuous self-similar solution gives no information for the stability of the actual 
dynamical solutions. 



V. SUMMARY 



We have studied the collapse of the Brans-Dicke scalar field to examine the dependence of the critical phenomena 
on the theory of gravitation. First, by conformally transforming the Roberts's solution into the Brans-Dicke frame, we 
find for lo > —3/2 that a continuous self-similarity continues and that the critical exponent does depend on lo. Second, 
by conformally transforming the Choptuik's solution into the Brans-Dicke frame, we find for lo > —3/2 that at the 
critical solution shows discrete self-similarity. However, the critical exponent depens strongly on lo, while the scaling 
parameter depends weakly on lo. This weak dependence on the theory of gravitation may justify the approximation 
employed by Price and Pullin. || They find that the critical solution can be well approximated by a flat spacetime 
scalar field solution. Numerical dynamical solutions show discrete self-similarity. These are another examples of 
discrete self-similar critical solutions. 

Given these numerical examples, the critical phenomena are not universal in the original sense. The apparent 
analogy to statistical system may be just analogy not fundamental. Rather, the phenomena numerically observed can 
be regarded as a manifestation of the features of the infinite dimensional dynamical system. From this point of view, 
there are many interesting problems to be solved. For example, "Does the discrete self-similarity corresponds to the 
limit cycle?" or, more generally, "How can we classify the gravitational collapse phenomena from the point of the 
dynamical system?" , "Are there any examples which give rise to bifurcation phenomena?" 

We hope to report on the resolution of these problems in the future. 
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FIGURE CAPTIONS 



Fig.l. The scalar field on the center with a marginally subcritical initial parameter as a function of logarithm of the 
central proper time r = — log(t* — t), where t* is the time when an infinitesimal black hole forms. In this time 
scale the scalar field oscillates with a periodicity A ~ 3.43. 

Fig. 2. The mass of a formed black hole for supercritical solutions as a function of p — p* with p* being the critical 
initial parameter. The least square fit gives the exponent (3 ~ 0.38. 

Fig. 3. The scalar field on the center for 2uu + 3 = 0.5, 1, 2 with a marginally subcritical initial parameter as a function 
of logarithm of the central proper time t = — log(t* — t). 

Fig. 4. The mass of a formed black hole for supercritical solutions for 2uj + 3 = 2 as a function of p — p* . The least 
square fit gives the exponent /3 ~ 0.43. 
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